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Spin qubits involving one or two spins have emerged as potential building blocks for quantum information processing 
applications, resulting in many double quantum dot (DQD) studies [lH4(. Coherent control of a two-electron spin 
qubit close to the singlet /triplet (S/T+) anticrossing through Landau- Zener- St iickelberg (LZS) oscillations [El, 0] has 
been studied theoretically [ 7H9 | and demonstrated experimentally [10[ in DQDs. Recent advances with triple quantum 
dot (TQD) technology [lll-[l7( have suggested additional advantages, such as their potential for encoding quantum 
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FIG. 1: Device, stability diagram, and rise time dependence of LZS oscillations, (a) Electron micrograph of a device 
identical to the one measured. Gates 1 and 2 are connected to high frequency lines for the application of fast voltage pulses 
(5Vi,5V2) in addition to DC voltages (Vi,V2). Gate C tunes the (1,1,1) region size, (b) Stability diagram obtained from 
numerically differentiating the left QPC detector conductance with respect to V2 at B=0.2 T. Black is low, red is medium, 
and yellow is high. Charge addition lines appear black, and charge transfer lines appear yellow [151 ]. A possible detuning line 
is drawn, (c) Schematic energy diagram illustrating the interplay resulting from the proximity of four of the three-electron 
spin states at finite B, including the hyperfine splitting, (d) Calculated pulse shapes for pulse duration r=10 ns after Gaussian 
filtering, leading to rise times of 6.6, 3.5, and 0.4 ns. (e-g) Left panel: experimental maps in the r-Vi plane from the numerical 
derivative of the right QPC conductance with respect to detuning along Vi for the pulse rise times corresponding to the 
colour code in (d). The pulse involves (5Vi,(5V^)=(4.0 r 1.7) mV to traverse the charge transfer line between (1,0,2) and (1,1,1) 
and repeats every r m =5 fis. B=45 mT. Black is low, orange is medium, and yellow is high transconductance. Right panel: 
calculated P A / maps in the r-Vi plane including the coupling to the (Ai/2,Qi/ 2 ) states for the same rise times as in the 
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left panel. Black (white) is low (high). Green [blue] dashed lines (guides to the eye) with negative [positive] slopes point to 
features pertaining to the Ai/ 2 /Q 3 / 2 [^1/2/^1/2^1/2)] interaction schematically represented in (c). 
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information [18j, that may soon be possible. Towards these goals we demonstrate, for for first time, the coherent 
manipulation of three-particle spin states in a TQD where all three spins play a role. 

The TQD is shown in Fig. QJa). Charge detection measurements are made with either one of the two quantum 
point contacts (QPCs) [19[ in the gate 1 and 2 voltage plane. High frequency pulses are applied via bias-tees (see also 
Supplementary Information Section I). 

We tune the device to the spin qubit regime using gate C [Fig. Ufa)] gaining control of the (Nl,Nc,Nr) = ( 1,1,1) 
electronic configuration between (1,0,2) and (2,0,1), where L, C, and R refer to the left, centre, and right quantum 
dots, respectively. The charge detection stability diagram obtained in the absence of pulses is shown in Fig. [ljb). 
We focus on coherent spin manipulation between the three-electron spin states of the (1,0,2), (1,1,1), and (2,0,1) 
electronic charge configurations. The dashed line illustrates a possible detuning line e, and e + [e_] corresponds to the 
location, along the detuning line, of the (2,0,1)/(1,1,1) [(1,0,2)/(1,1,1)] charge transfer line. 




FIG. 2: Three-electron spin state energy diagrams and spin funnels and arches, (a and c) Calculated energies 
vs. detuning e for two different sizes of the (1,1,1) region, where three spin-1/2 electrons occupy the triple quantum dot in the 
presence of a Zeeman splitting E z . The detuning line is at a 45° angle in the V1-V2 plane, (b) Left [right] panel: Spin funnel 
(half of the spin arch) in the numerical derivative of the right QPC conductance with respect to Vi [left QPC conductance 
(with a plane subtracted)] mapped in the detuning-B plane for the case of a wide (1,1,1) region. (<5Vi,<5V2)=(4.0,-1.7) mV 
[(-3.75,6.6) mV] traversing the charge transfer line between (1,0,2) [(2,0,1)] and (1,1,1). The pulse has a t=100 ns duration 
(>T 2 ) and repeats every r m =5 {is. The DC detuning axis is purely along Vi [V2]. (d) Numerical derivative of the left QPC 
conductance with respect to V2 in the presence of a pulse with (£Vi,£V2)=(-8,8) mV across the charge transfer line between 
(2,0,1) and (1,1,1) for the case of a narrow (1,1,1) region. The spin arch is fully seen in this case. r=50 ns and r m =2 fis. The 
DC detuning line is at 45° in the V1-V2 plane, permitting both sides of the spin arch to be observed. The dashed lines in (b) 
and (d) are theoretical fits [with detuning-dependent interdot couplings in Fig. [2jb)]. The experimental e axes in (b) and (d) 
are offset by the appropriate pulse length with respect to (a) and (c). 

We focus our attention on four of the three-electron spin states, as shown in Fig. Die). One of the lowest energy 
eigenstates of the system is a spin doublet with total spin 5^=1/2 and S z =+l/2, labelled A^ 2 following Ref. jl7j . 
The two lowest energy quadruplet states with S z =+%/2 and +1/2 are labelled Q 3 / 2 and Qx/2, respectively. A higher 
energy doublet state, A 1/2? lies m close proximity to Qi/ 2 - 

Filtering a 10 ns rectangular pulse controls the rise time [Fig. QJd)]. At long (short) rise times, pulses appear Gaus- 
sian (almost rectangular). When the detuning pulse is applied such that the system passes through the A / 1 ^ 2 /Q 3 / 2 an- 
ticrossing, mixing between all four states can occur due to a non-zero Landau-Zener transition probability (Supple- 
mentary Information Section VI). The speed through the anticrossing (related to the rise time) determines the final 
superposition of quantum states. 

Applying a detuning pulse of duration r will result in a phase accumulation between the quantum state components. 
This phase is related both to the accumulation time and the detuning voltage (El-Hof. The resulting LZS oscillations 
in the A^ 2 state final projection probability are measured using standard spin-to-charge conversion techniques (20j . 

The dependence of the LZS oscillations on rise time is shown in the left panel of Fig. DJe-g) for a large enough 
detuning pulse to allow mixing with the (Ai/ 2 ,Qi/2) states. At long rise times [Fig. QJe)], LZS oscillations from the 
A^y 2 /Q 3 / 2 qubit are observed. The maxima in the signal have a negative slope in the r-Vi plane (see green dashed 
lines for a guide to the eye). The LZS oscillations disappear progressively as the rise time decreases, since the Landau- 
Zener tunnelling probability approaches unity. A new set of oscillations emerges at very short rise times [Fig. Hg)]. 
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These oscillations have a positive slope and a long period (blue dashed lines). Calculations confirm that these long 
period oscillations are due to the A / 1 ^ 2 /(A 1 / 2 ,Qi/2) interactions [see Fig. [TJc)]. The calculated results including the 
additional interactions (see Supplementary Information Section III) agree well with experimental data and are shown 
in the greyscales of Fig.QJe-g). For more details about the interplay between coherent processes involving the different 
states as well as the Bloch sphere repres 

The calculated energy diagram for all the spin states of the three-electron system in the absence of hyperfine coupling 
but in the presence of Zeeman splitting E z from an in-plane magnetic field B is plotted versus detuning in Fig. [2fa,c) 
oilOl for (1,1,1) regions of different sizes with different couplings. The size of the (1,1,1) region is readily tuneable 
|15j . The quadruplets (with 5^=3/2) are defined to be independent of detuning. The high energy antibonding 
doublet states are also shown in grey in Fig. EJa,c). More details about the three-electron spin Hamiltonian can be 
found in Supplementary Information Section II and in f\j\ . 

The exact location of the hyperfine-induced anticrossings between A^ 2 and Q 3 / 2 along the detuning axis [open 
circles in Fig. EJa)] depends upon the Zeeman splitting, which results in a B dependence. This is shown for the two 
A^/ 2 /Q 3 /2 anticrossings in Fig.[2fb), where the individual spin funnels [l| measured along two detuning axes combine 
to form a "spin arch" across the expanded detuning range. For the case of the narrow (1,1,1) region, the spin arch can 
be measured in a single detuning line running across (1,1,1) and detected by a single charge detector. This enables 
coherent manipulation between three-electron spin states, as the spin arch contains the exchange energy contributions 
from both relevant combinations of right-centre and left-centre dots. Fits to the model described in Supplementary 
Information Section III are superimposed with the data in Fig. |2fb,d). 
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FIG. 3: LZS oscillations from the two A1/2/Q3/2 q u °its for a wide (1,1,1) region, (a) [(d)] Numerical derivative 
of the right [left] QPC conductance with respect to detuning along Vi [V2] illustrating oscillations vs pulse duration r across 
the (1,0,2) [(2,0,1)] to (1,1,1) charge transfer line at B=60 mT, (£Vi,£y 2 ) = (4.0,-1.7) mV [(-4.11,7) mV]. In (a), both V 2 and 
Vi are swept in order to detune parallel to the pulse in the V1-V2 plane. (b,c and e,f) Probability of ending in the A^ , 2 state 

as a function of r taken at two consecutive detunings further away from the crossing and fit for T 2 . In (b,c), the pulse goes 
from (1,0,2) to (1,1,1) while in (e,f), the pulse goes from (2,0,1) to (1,1,1). The experimental data [fit] are shown as points 
[theoretical lines colour-coded by the filled circles in Fig. [Ha)]. The values of T 2 extracted from the single parameter fit to the 
LZS model are indicated. 

In the case of the large (1,1,1) region, for r<25 ns we observe coherent LZS oscillations as a function of detuning and 

r, as shown in Fig.[3ja,d). The probability of measuring the A^ 2 state after the pulse, P^/ , is plotted for two pulse 

1/2 

detunings. From a single parameter fit (see Supplementary Information Section III), we extract values for T 2 , while 
the position of the fringes is completely determined by independently measured parameters (not shown) following 
methods from Ref. [10]. The frequency of the LZS oscillations is directly related to the level spacing. This manifests 
itself in the data with a decreased frequency the closer the pulse is to the anticrossing, as expected. In general (not 
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shown) such fits reveal that longer T 2 values are obtained for pulses reaching further past the anticrossing. We believe 
such a behaviour would be consistent with charge noise effects due to the flattening of A^ 2 vs. detuning. Data for 
the magnetic field dependence of the LZS oscillations are in Supplementary Information Section V. 




1 6 11 16 1 6 11 16 

X (ns) x (ns) 



FIG. 4: Coherent three-electron spin state manipulation for a small (1,1,1) region, (a) Stability diagram for a small 

(1,1,1) region, obtained from numerically differentiating the left QPC detector conductance with respect to V2 at B=60 mT. 

(£Vi,£V2) = (-8,10) mV to traverse the charge transfer line between (2,0,1) and (1,1,1) and reach the (1,0,2) region. r m =2 {is. 

Black is low, orange is medium, and yellow is high, (b) Higher resolution data zooming into the spin-to-charge conversion 

region (2,0,1) to show the details of the coherent features observed as parallel lines, (c) Experimental map in the T-V2 plane 

revealing the coherent behaviour of the three-electron spin states from the numerical derivative of the left QPC conductance 

with respect to V2. Vi is swept proportionally to V2 in order to detune parallel to the pulse direction. The spin-to-charge 

conversion is performed in the (2,0,1) region for all the coherent oscillations observed in (c). (d) Calculated Pa/ maps in the 

^1/2 

r-e plane for the same experimental settings as in (c). 

In the case of a narrow (1,1,1) region, as seen from the stability diagram in Fig. Ufa), a set of oscillations parallel to 

the opposite charge transfer line (1,0,2)-(1,1,1) are observed in the (2,0,1) spin-to-charge region when the pulse is large 

enough to engage the third spin in the coherent modulation by reaching the second (i.e. far) A^ 2 /Q 3 / 2 anticrossing. 

These new oscillations involve the coherent spin-state manipulation of three spins simultaneously. The zoomed-in 

version of these oscillations in the stability diagram is shown in Fig. Hfb). With a detuning line parallel to the pulse 

direction in the stability diagram (i.e. parallel to the center dot addition line), it is possible to map out the coherent 

behaviour among the three-electron spin states, as seen in Fig.|4jc). The set of closely spaced oscillations with negative 

slopes in the left part of the r-V 2 plane corresponds to the LZS oscillations from the A / 1 ^ 2 /Q 3 / 2 qubit that is close 

to the (2,0,1) region, while the narrow curved features in the right hand side of the r-V 2 plane correspond to LZS 

oscillations involving the A^ 2 /Q 3 / 2 and A / 1 ^ 2 /(A 1 / 2 ,Q 1 / 2 ) interactions close to the (1,0,2) region. The corresponding 

theoretical P^y map is shown in Fig. 21(d). The two broad fringes marked by arrows are due to the A^ 2 / (Ai/ 2 ,Q]y 2 ) 
1/2 

interactions, while the well-defined fringe above them is due to the A^ / 2 /Q 3 / 2 interaction. Above this is a dense group 

of fringes, where the relevant energy splittings grow sharply towards the (1,0,2) region producing very fast oscillations 

that are resolution limited by pixelation. See Supplementary Information Section VIII for an example of measured 

and calculated Pa/ maps in the (pulse amplitude)-detuning plane showing the relevant pulse amplitudes where 
^1/2 

coherent rotations from the two A' x / 2 /Q 3 / 2 qubits are observed. 

In conclusion, we have performed coherent spin manipulation in the form of Landau-Zener-Stiickelberg oscillations 
for both A / 1 ^ 2 /Q 3 / 2 qubits in a triple quantum dot. The observation of the spin arch and the resolution of coherent 

* 

oscillations from the three-electron spin states compare well with calculation results. Fits permit the estimation of T 2 , 
which lies in the 5-18 ns range (18 ns not shown). The coherent interplay between A^ 2 /Q 3 / 2 and A^ 2 / (Ai/ 2 ,Q]y 2 ) is 
required to understand the results. The ability to coherently mix multiple three-spin states combined with micro- 
magnet technology is expected to lead to the universal control of such states. Further improvements of coherent spin 
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manipulation techniques in triple quantum dots will allow more complex pulse sequences to achieve arbitrary coherent 
rotations in the Bloch sphere of a qubit encoded in three-electron spin states. 
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I. METHODS 



The device is fabricated on a GaAs/AlGaAs heterostructure grown by molecular beam epitaxy with a density of 
2.1 xlO 11 cm -2 and a mobility of 1.72xl0 6 cm 2 /Vs. Ohmic contacts are used to contact the two-dimensional electron 
gas (2DEG) located 110 nm below the surface. TiAu gate electrodes are patterned by electron-beam lithography to 
allow electrostatic control of the triple quantum dot (TQD). Two gates are used to define quantum point contacts 
(QPCs) used as charge detectors on the left and right of the TQD. 

Charge detection measurements are made by measuring either the left or right QPC conductance with a lock-in 
technique using a typical root-mean-square modulation in the 0.05-0.1 mV range. The QPC detector conductance 
is tuned to below 0.1 e 2 /h. High frequency pulses from two synchronized Tektronix AWG710B are applied via a 
bias-tee. The pulse of duration r is typically <25 ns and the waveform is typically repeated every 2 to 5 /is. The 
pulse rise times are controlled by passing the programmed rectangular pulses through low-pass filters internal to the 
AWG10B or through external Mini- Circuits SBLP filters. Typical rise times are 6.6 ns. The device is bias-cooled 
with 0.25 V on all gates. Once cold, suitable gate voltages are applied to the gates to form the TQD potential. The 
dilution refrigerator has an electron temperature of ~110 mK, as determined from the temperature dependence of 
the lineshape for an addition line. 
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II. THREE-ELECTRON SPIN HAMILTONIAN 

Based on Ref. [17j, the Hamiltonian for a system of three electron spins in the presence of a magnetic field along z 



is: 



H = J LC j Sl • Sc — - 



Jrc [ Sr - Sc — - 



E z {Si 



(1) 



where is the exchange interaction between spins in dots i and j, Si is the spin in dot i and ^ is the Zeeman 
energy. According to Ref. [17j, the three-spin system is characterized by eight eigenvectors, which are divided into 
two subgroups by the exchange energy: four quadruplet states Q with a total spin S = 3/2 (S z = ±3/2, ±1/2) and 
two pairs of doublet states A and A' with a total spin S = 1/2 (S z — ±1/2). We refer the reader to Ref. [17j for the 
complete list of eigenstates and eigenvalues, and we write down only those that played a role in the main text: 



\q 1/2 ) = -L(itu) + mt> + litt)) 



i 



|Ai /2 > = ; 

1 V 4 " 2 + 2ft ( j lc ~ 2Jrc) 

IQ3/2) = mt> 

K/ 2 > = 1 

1 V4^ 2 + 2^(2 J RC - Jlc) 



{{Jlc - Jrc + fi)|m> + (Jrc - fi)|tit> - Jlc\IK)) 



((-Jlc + Jrc + fi)|tU> - (Jrc + mm) + J L c\itt)) 



(2) 

(3) 
(4) 
(5) 



where ft = yj J\ c ± J\ c — Jlc Jrc and the eigenvalues are: 



^1/2 
E A 1/2 



E, 



3/2 



l l/2 



-E z /2 

-(Jlc + Jrc ~n + E z )/2 
SE z /2 

-(Jlc + Jrc + £1 + E z )/2 



(6) 
(7) 
(8) 
(9) 



In the limits of large |e|, the A^ 2 (A^) doublet state evolves to the D^ 2 (D 1 / 2 ) and D\j2 (^1/2) states of 
Ref. [13 • 

III. MODEL FOR LANDAU-ZENER-STUCKELBERG OSCILLATIONS IN THE A' 1/2 /Q 3/2 SYSTEM 

Following the terminology of Ref. [17j the Hamiltonian for the Landau-Zener-Stiickelberg (LZS) model in the 
A i/2/Q3/2 system is: 



H = 



r* A , 



3/2 



r A',Q 3 



A / ,Q 3 /2 



(10) 



x l/2 



where the off-diagonal term Fa / ,q 3/2 is the A / 1 ^ 2 /Q 3 y 2 coupling originating from the hyperfine interaction betwen 
the electron spins and the nuclear spins via the x and y components of the Overhauser field gradients between the 
dots. This Hamiltonian is equivalent in the limit of weak "spectator dot" coupling to the 2-dot Hamiltonian in the 
S/T + basis described in Ref. 

For situations in which the (Q 1 / 2 ,Ai/2) states play a role we use a Hamiltonian of the form: 



/ e q r A , Ql 



H = 



1/2 

r A,Q 1/2 E A 1 



/2 



\ X A',Q 1/2 



/2 


1 A' A 



Q3/2 
1 A',Q 3/2 



r A',Q 1/2 \ 

r A',Q 3 /2 
^1/2 / 



(11) 
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The (rA',Q 1/2 ,rA',A?rA,Q 1/2 ) couplings are due to z component of the Overhauser field gradients between the dots. 
This Hamiltonian is equivalent in the limit of weak "spectator dot" coupling to the 2-dot Hamiltonian in the S/T+/T0 
basis. 

Note that the magnitudes of the off-diagonal coupling elements are empirically fitted to the observed magnitude of 
the LZS oscillations. This has no significant effect on the period of the LZS oscillations. The couplings are typically 
- 0.1-0.2 peV. 

The time evolution of the density matrix p is calculated from the initial state at large detuning where probability 

Pa/ =1, as described by the following equation: 
^1/2 

< f t =i[p,H/h] (12) 

The solution of the time evolution of p involves a series of differential equations solved numerically by the Runge- 
Kutta method. To simulate decoherence effects appropriate off-diagonal terms are included in the derivative of 
the density matrix leading to exponential decay of the resulting oscillations. The pulse shape is simulated by the 
convolution of a rectangular pulse of length r with a Gaussian -j=-e~ l / 2s where s is the Gaussian time, which is 

approximately equal to half of the measured rise time from 10% to 90%. The applied magnetic field and the nuclear 
field gradient from the difference in Overhauser fields are kept constant. At the end of the pulse the density matrix 
in the (Q1/2 5^1/2 5Q3/2 ^1/2) basis is projected back onto A^ 2 corresponding to the charge measurement. 



IV. STABILITY DIAGRAMS WITH LONG r PULSES 



For a pulse with long duration r=100 ns, a line, corresponding to the measurement location where the end of the 
pulse reaches the A^ 2 /Q 3 / 2 anticrossing, appears inside the stability diagram (black triangle), not far from the charge 
transfer line (black circle) inside the respective spin-blockade regions of (1,0,2) and (2,0,1) (Fig. [5]). We call this new 
line the Ay 2 /Q 3 / 2 line. In the absence of a pulse, we observe only the charge transfer line, i.e. the boundary between 
the (2,0,1) and (1,1,1) regions. 




FIG. 5: (a) Numerical derivative of the right QPC conductance with respect to Vi in the presence of a pulse with 
(£Vi,£V2) = (4.0,-1.7) mV across the charge transfer line between (1,0,2) and (1,1,1). Black is low, red is medium, and yel- 
low is high transconductance. The pulse has duration t=100 ns (>T 2 ) and repeats every r m =5 ps [the repetition rate is chosen 
< Ti, which is ~10 ps in this system (not shown)]. B=60 mT.(b) Numerical derivative of the left QPC conductance with 
respect to V2 in the presence of a pulse with (£Vi,£V2) = (-3. 75,6.6) mV across the charge transfer line between (2,0,1) and 
(1,1,1). r=100 ns, r m =5 ps, and B=83 mT. 



V. MAGNETIC FIELD DEPENDENCE 



Figure [6] contains the investigation of the magnetic field dependence of the LZS oscillations. The number of 
oscillations grows with B, as the energy difference between the two states in the qubit grows with E z due to the 



8 



curvature in A' x , 2 . These experimental results compare very well to the calculations made with the single qubit 
model also shown in Fig. [6] This completes the characterization of the A / 1//2 /Q 3 / 2 qubit. The truncation of the LZS 
oscillations along the Vi axis in the experimental data is due to a small spin-to-charge conversion region in this case. 




FIG. 6: Magnetic field dependence of the LZS oscillations. White (black) is low (high) transconductance. From left to right, 
B[mT] = 10, 30, 45, 60, and 90. The top row contains the experimental maps in the r-detuning plane for the numerical derivative 
of the right QPC conductance with respect to detuning along Vi . The bottom row contains the maps calculated from the LZS 
model of Eqn.HJB White (black) is high (low) Pw . 



VI. LANDAU-ZENER PROBABILITY 



In order to measure the transition probability Plz for the avoided crossing between states A^ 2 and Q 3 / 2 , we 
adopt the following procedure, originally used in a DQD system (lp) . From the (2,0,1) configuration, we induce the 
(1,1,1) transition by going non-adiabatically through the avoided crossing between the A[ , 2 and Q 3 / 2 in 0.5 ns. In 
this way, the A[ , 2 is preserved during the sweep. Immediately after that, we sweep through the avoided crossing in 
the opposite direction and measure the probability of returning in A^ 2 , P^ , as a function of the return sweep 

time. Figure [3 shows the results. When the sweep time for the return is short, i.e. in the non-adiabatic regime where 
AE/v — >• 0, Landau-Zener tunneling is efficient so Plz = 1 (see Ref. [10]). As the sweep time is increased, Plz 
shrinks so the probability of finding the system in Q 3/2 increases, which appears as an exponential decay of Pa/ . 

1 ^1/2 

The characteristic time obtained by fitting an exponentially decaying function to the data is of 260 ns. This implies 
that to obtain a superposition of equal weights between A! x , 2 and Q 3 / 2 during a single sweep through the anticrossing, 
it is necessary to sweep through the transition for close to but less than 260 ns. This is impractical, as the decoherence 
time T 2 is close to 15 ns. 



1.0 




Sweep time (jlis) 

FIG. 7: Landau-Zener transition probability for the avoided crossing between A[/ 2 and Q 3 / 2 as a function of sweep time. The 
characteristic time extracted from the exponential fit is 260 ns. 
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VII. INCLUDING HIGHER ENERGY LEVELS AND BLOCH SPHERE TRAJECTORIES 




M |a«) K) 



FIG. 8: Probability P A / measured at fixed detuning. Colored solid lines represent the best fit using (a) the single qubit 

^1/2 

model and (b) including A / 1 / 2 /(A 1 / 2 ,Q 1 / 2 ) interactions, (c) through (e): qubit trajectory on the Bloch sphere for the first 
three extrema in (a), i.e. for 7r, 27t, and 3ty accumulated phase. The red part of the trajectory starting at A[, 2 occurs before 
the pulse reaches the anticrossing. The dark blue curve corresponds to the rotation when the pulse is beyond the anticrossing 
(creating a coherent superposition), while the light blue portion is the trajectory after the pulse. No dephasing is included for 
the trajectories, so the trajectory remains on the surface of the Bloch sphere. All other parameters are as in (a). 



Including higher energy spin states in the model is required to fully describe the oscillations in P A / . The effect 

^1/2 

of the oscillations arising from the interactions to higher energy spin states is observed even with a 6.6 ns rise time, 
as demonstrated by the fits to the simple LZS model [Fig. Ufa)] and multi- level interacting model [Fig. Efb)]. Such 
an interplay is in accordance with theoretical predictions from Ref. 0|. 

It is possible to represent the qubit evolution on the Bloch sphere for the oscillations modelled in Fig. Ufa). We 
present three cases, where the accumulated phase during the pulse is 7r, 27r, and 3tt [see Fig. Efc), (d), and (e), 
respectively]. These trajectories represent the evolution of the qubit state vector for values of r corresponding to 
the first three extrema in Fig. 15(a). During the pulse, (dark blue portion of the trajectory), the projection of the 
state vector on the x-y plane makes 0.5, 1, and 1.5 rotations, respectively. Note that these trajectories do not have 
a superposition between AJ/ 2 and Q 3 / 2 below the equator of the Bloch sphere, as the pulse rise time for such a 

superposition would be much greater than T 2 [Tq| . 



VIII. MORE ON COHERENT SPIN MANIPULATION OF THREE-ELECTRON SPIN STATES 

It is possible to decouple the r dependence from the pulse amplitude changes that occur at small r (reduction of the 
rectangular pulse amplitude by Gaussian convolution), by fixing the value of r (e.g. 10 ns) and stepping the detuning 
pulse amplitude. In such a map, the location of the two A / 1 ^ 2 /Q 3 / 2 anticrossings are observed as two oblique yellow 
lines, seen in Fig. [9fa). The corresponding calculated map of P^y is shown in Fig. [9^b). 
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FIG. 9: (a) Experimental map in the (pulse amplitude)-V2 plane at r=10 ns revealing the coherent behaviour of the three- 
electron spin states from the numerical derivative of the left QPC conductance with respect to V2. Vi is swept proportionally 
to V2 in order to detune parallel to the pulse direction. The oblique yellow lines correspond to the two A^/Qg^ anticrossings. 
(b) Calculated P A / maps in the (pulse amplitudeW plane for the same experimental settings as in (a). 



